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ABSTRACT 

We use the distribution of extrasolar planets in circular orbits around stars with surface convective 
zones detected by ground based transit searches to constrain how efficiently tides raised by the planet 
are dissipated on the parent star. We parameterize this efficiency as a tidal quality factor (Q*). We 
conclude that the population of currently known planets is inconsistent with < 10^ at the 99% 
level. Previous studies show that values of between 10^ and 10^ are required in order to explain 
the orbital circularization of main sequence low mass binary stars in clusters, suggesting that different 
dissipation mechanisms might be acting in the two cases, most likely due to the very different tidal 
forcing frequencies relative to the stellar rotation frequency occurring for star-star versus planet-star 
systems. 

Subject headings: convectiov — planet-star interactions — stars: interiors — stars: rotation — stars: 
winds, outflows — turbulence 
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1. INTRODUCTION 

Exoplanets with orbital distances < 0.1 AU from 
their host stars, called close-in exoplanets, have pre- 
sented an especially puzzling challenge to theories of 
planet formation. The protoplanetary disk is too warm 
( 2000 K) so close to a star to allow the condensation 
and accumulation of icy an d rocky material required 
to form planets (ILin et all [1996. : \Miller et al.. .2009; 
llbgui fc Burrowsll2009| ). 

For close-in exoplanetary systems, their mutual tidal 
gravities induce significant tidal bulges in the planets 
and stars. Dissipation of the accompanying tidal energy 
drives obliquities to zero and rotation rates to near 
synchronous, processes that probably take millions of 
years for planets but billions of years for stars. While 
its orbit is eccentric, dissipation of tidal energy within 
a planet can reduce the orbital semi-major axis and 
eccentricity, as well as warming the planet's interior 
with significant co nsequ ences for th e planet's thermal 
evolution (Jackson et al.ll200 8b: Ibg ui fc Burrows! 120091 : 
iLiu et all 12008; .Miller et al.i ,2009 ). The majority of 
planet-hosting stars rotate more slowly than their 
close-in planets revolve, and so tides raised on these 
stars also reduce eccentricities and semi-major axes. 

Although the effects of tides raised on close-in ex- 
oplanets become negligible as eccentricities shrink, 
tides on host stars continue to reduce semi-major axes 
long after eccentricities are negligible, as long as the 



stellar rotation rate is smaller than the orbital mean 
motion. For the systems for which the rotation of 
the star is observationally constrained, the stellar spin 
period is known to be longer than the orbital period. 
Typically, planet-hosting stars are older than I Gyr, 
and so stars without reported rotation rates likely 
rotate slower than their planets revolve. Moreover, 
observational biases favor detection of planets around 
slow rotators. As a result, the tides on the host star 
tend to dominate the long-term tidal evolution of 
close-in planets. Eventually, the planets may cross 
their Roche limits (0.007 AU for a Jupiter-like planet 
around a Sun-like star), where they are tidally disrupted. 

On the other hand, tidal spin-up may synchronize 
the stellar rotation to the orbital period, in which case 
the planet will eventually reach a stable orbit. The 
total angular momentu m of the system determines 
which scenario occu r s (iCounselmanI 119731 : iGreenber^ 
I1974D . iLevrard et"all ()2009D show that all systems with 
transiting planets found to date by ground transit search 
surveys except HAT-P 2 b have insufficient angular 
momentum to prevent this destruction. Loss of angular 
momentum through shedding of stellar wind dooms 
even HAT-P 2 b. Thus, given enough time, the loss 
of close-in exoplanets through orbital decay is inevitable. 

Tidal evolution of an orbit increases rapidly for 
decreasing semi-major axis, and so the distribution of 
semi-major axes for observed planets is sensitive to 
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the rate of tidal dissipation. Since the probability for 
a planet to transit its host star increases for planets 
nearer their stars, transiting planets are especially 
susceptible to tidal effects. The rate of orbital decay 
and frequency of tidal destruction also depends sensi- 
tively on the rate of tidal dissipation within the host 
star. This rate is r elated to the efhciency parameter 
(jGoldreichl I1963D : larger corresponds to less 
tidal dissipation and slower orbital evolution. The 
origins of tidal dissipation within gaseous planets and 
stars have been studied extensively, but remain poorly 
understood, with estimates based on theoretical and 
observation a l stu d ies ranging from 10^ to 10 ^ (c.f . 
ZahnI [19661 fl970l ITOTI IiQTTI: IZahn fc Bouchetl [1989 : 
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Consequently, the time until a planet crosses its Roche 
limit and is removed, which we will call a planet's time 
left (TL), depends both on and on its current semi- 
major axis, among other parameters. For a population 
of tidally evolving planets, we expect to find few planets 
with TL much less than the whole lifetime or current 
system age. Otherwise, we would conclude that we 
have caught a large fraction of the planets in the last 
extremely short moments of their lives, just before they 
are disrupted by their star. By tuning until we 
generate a statistically likely distribution of TL, we can 
constrain and the frequency of tidal disruption of 
exoplanets. 

As discussed in Section 2, observational biases have 
important and complex infiuences on the distribution 
of calculated TL- values for transiting planets and must 
be considered in order to produce statistically reli- 
able co nstraint s on Q„. Several previous studies (e.g., 
iCarone fc Pa tzoldl 120071 ) have attempted to place con- 
straints on using considerations similar to ours but 
only applied to individual planets. Results from some 
of those studies are consistent with > 10'', but in- 
ferences based on the orbital evolution of a single planet 
may not be statistically meaningful. 

The outline of this paper is as follows: in Section [2] 
we describe our methods and assumptions for calculat- 
ing the orbital evolution of planetary systems; in Section 
[3] we show the sample of transiting planets and the cor- 
responding parameters which were used in this work; in 
Section 2] we discuss the various observational and as- 
trophysical biases that affect the sample of systems with 
transiting planets found by ground based transit searches 
and our procedures and assumptions for how to deal with 
them in the analysis; in Section [5] we outline the proce- 
dure we use to derive constraints on the value; in 
Section|n]we show our limits to Q,; in Section[7]we show 
the apparent discrepancy between our results and esti- 
mates of derived from binary stars in open clusters. 



2. ORBITAL EVOLUTION 



Tidal decay of close-in planets involves the exchange 
of angular momentum between a planet's orbit and its 
host star's rotation. For the stars we consider here, sev- 
eral processes influence the stellar rotation, in addition 
to tidal processes, and accurate modeling of the orbital 
decay requires consideration of these effects. For this 
purpose, we solve the following system of ordinary dif- 
ferential equations: 
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where M* is the mass of the star; i?, is the radius 
of the star; nip is the mass of the planet; is the 
tidal quality factor of the star; sign(LLiconv — i^orb) takes 
the value 1 when the stellar convective zone is spinning 
faster than the planet and -1 when the reverse is true; 
K — 0.35 M0i?Qday^Gyr~"'^ is the proportionality con- 
stant, parametrizing the strength of the magnetic wind 
of the star; Wgat = 1-84 day~^ is the wind saturation 
frequency; /conv is the moment of inertia of the stellar 
convective zone; Lconv is the angular momentum of the 
stellar convective zone; /rad is the moment of inertia of 
the stellar radiative core; Lrad is the angular momentum 
of the stellar radiative core; Tc = 5 Myr is the stellar 
core-envelope coupling timescale; Mj-ad is the mass of the 
stellar radiative core; i?rad is the radius of the radiative- 
convective boundary in the star; Wconv = Lcom/ Icom is 
the angular frequency of the stellar convective zone. 

We wish to follow a planet-star system as its semi- 
major axis shrinks under the influence of tidal friction. 
There are two sources of friction: the tides on the star 
and those on the planet. However, the latter is only im- 
portant as long as either the orbit is eccentric and/or 
the planet is rotating asynchronously. Since the angu- 
lar momentum stored in the rotation of the planet is 
quite small compared to the orbital or stellar spin angu- 
lar momenta, it is safe to assume that the planet spin 
is synchronized quickly compared to any orbital evolu- 
tion. Further, we will restrict our sample to only systems 
with nearly circular orbits. In this case, the evolution of 
the semi-ruajor axis (a) is given by Equation fl]) above 
(|Goldreichlll9"6aiKaulalll968tlJackson et al.ll2008ari . This 
expression makes the approximation that the planet's 
mass can be neglected compared to the star's mass, a 
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perfectly reasonable assumption for all the systems we 
consider, given the uncertainty in the value of Q*. 

The angular momentum that is taken away from the 
orbit by the tidal friction is deposited in the star, acting 
to spin it up, and while for most of the known tran- 
siting planets the orbit does not have enough angular 
momentum to spin up the star to synchronous rotation, 
for at least one planet this is not true. Further, since 
the tidal friction couples the planet with the surface 
convective zone of the star, it is possible that if the 
core-envelope coupling is not sufficient, synchronous 
rotation can be imposed on the envelope only. For 
this reason we also follow the evolution of the angular 
momentum of the star {L^,) - Equation ([2]). Here, we 
do not make the approximation M^, ^ mp, like we did 
for the orbital evolution above, in order to make the fi- 
nal set of equations conserve angular momentum exactly. 

In addition, since we are interested in timescales of 
order Gyrs, we cannot neglect the spin-down of the 
host star due to its own magnetic wind. The effects of 
stellar wind shedding on rotation have been extensively 
studied, but remain poorly understood. Equation ([5]) 
represents the current best description of these effects 
and is mot ivated by a combination of theory and 
observation (iStauffer fc Hartmannlll987t iKawalerl [19881: 
IBarnes fc Sofialll996D . Note that we do not introduce 
any free parameters to describe the stellar wind. The 
value of K in Equation ([5]) is determined by the present 
rotation rate and age of the Sun and the value of the 
wind saturation frequency Wsat comes from fitting the 
observations of stellar rotation rates in open clusters of 
different ages. With those values, the expression above 
matches the rotational evolution of single stars in open 
clusters at young ages, as well as the rotation rate of 
the Sun and other less well constrained but older stellar 
populations. Hence, theoretical interpretations aside, it 
can be viewed as a parameterization of the rotation rate 
observations over the range of ages we encounter during 
the orbital evolution of the exoplanet systems in our 
sample. 

Open cluster rotation rates impose one more com- 
plication on our model: core-envelope de-coupling. 
In order to explain how stars with a wide spread 
of rotation rates at young ages end up with similar 
rotation rates later, it is necessary to allow for quickly 
rotating cores in the slow surface rotation rate stars, 
and a re-distribution of the excess angular momentum 
to the convective zone at a later time (llrwin et al.l 
[2007t llrwin fc BouvIctI [2009t iDenissenkovl I2010D . The 
expressions for th e separ a te cor e and envelope evolution 
were derived by lAllairj (119981 ). and for planet hosts 
take the form of Equations ^ and ([6]). We test the 
sensitivity of our results on the assumed core-envelope 
coupling by repeating our analysis under the assumption 
that the core and envelope are perfectly coupled (tc = 0). 

Finally, Equation ([7]) for the st ellar core-envelop e 
differential rotation was proposed bv lMacGregod (|1991l ). 

Note that in Equation ([ij we treat as having a 
fixed magnitude. We feel that more complex assump- 
tions are not justified, since the dependence of on 



frequency is poorly known, and there is substantial 
disagreement between observational and theoretical 
estimates as discussed in the introduction. 

In this work we start a planetary system's evolution 
from an age of 5 Myr after the (model) birth of the 
star, with an initial orbital separation that evolves to 
the observed semi-major axis at the present system age. 
The initial age is assumed to be early enough so that no 
significant spin up of the star due to the tides raised by 
the planet has occurred, but any non-tidal evolution of 
the planetary orbit has stopped. This is reasonable if 
one assumes that hot Jupiters arrive at their extremely 
close-in orbits through disk migration, since by that time 
the p rotoplan etary disk ha s dissipated ([Haisch et al.l 
120011 : [Bouwnia n et al.l 120061 ). but for other migration 
mechanisms this might not be the case. 

The advantage of starting the evolution at 5 Myr is 
that with these assumptions the initial stellar rotation 
distribution is relatively well constrained from obs erva- 
tions of young open clusters (jlrwin fc Bouviedl200l and 
references therein). However, this age is significantly 
before the star has arrived at the main sequence, so all 
stellar parameters, in particular the radius, the masses 
of the radiative core and convective envelope, and 
the corresponding moments of inertia all significantly 
evolve with time. This forces us to allow for stellar 
evolution along with the orbital and spin evolution of 
the planet-star systems. In order to follow the evolution 
of the star, we use evolution tra cks calculated using the 
YREC stellar evolution code (jPemargue et al.l I2008D 
with masses OAMq, 0.5Mq, 0.6Mq, O.TMq, O.8M0, 
0.9Mo, I.OMq, 1.O5M0, l.lMg and 1.25Mo. The 
evolution of individual stars is determined from cubic 
spline interpolation within this grid of models. Stars 
below OAMq and above 1.25Mo are excluded from our 
analysis (see below). 

In our evolution we assume that the tides couple the 
orbit of the planet only to the convective zone. While, 
in general, we should split the tidal torque in two parts, 
one spinning up the convective zone and the other 
spinning up the core, the coupling to the core is likely 
to be negligible compared to the convective zone for 
two reasons: (1) the amplitude of the tidal deformation 
scales strongly with radius, so it will be much smaller for 
the core than for the convective zone; (2) currently there 
is no known mechanism for dissipating the tidal energy 
in the core, so averaged over an orbit there should be 
no net angular momentum transfer from the tides to the 
core. 

An additional complication for our models is that 
the prescription described above is only valid for low 
stellar masses. For masses larger than approximately 
1.2Mq the surface convective zone becomes negligible 
in mass, so we cannot treat it as the only s i nk for 
angular momentum. Further, iWolff fc SimonI (|1997l ) 
indicate that the angular momentum loss described 
above is only valid for low-mass stars (M* < 1.3Mq). 
For more massive stars, presumably the lack of surface 
convection suppresses the stellar wind, and the loss of 
angular momentum is much weaker. For reasons of 
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numerical stability, for stars above I.IMq we ignore 
the core-envelope decoupling, treating the star as a 
solid body and we completely exclude from our analysis 
stars above 1.25Mq, since those stars do not posses a 
significant surface convection zone and could be subject 
to a different mode of tidal dissipation, in addition to 
weaker stellar wind. Consequently, close-in planets 
around very massive stars may have much larger TL 
values than around less massive stars, with possible 
implications for planet surveys of massive stars. 

One particular effect of including the stellar wind 
for low-mass stars is that even after the star has 
synchronized its spin with the orbit, the semi-major 
axis continues to decay. In fact, a positive feedback 
loop is created at this point: tides keep the stellar spin 
synchronized with the planet, so as the stellar wind 
removes angular momentum, the orbit continues to 
shrink. This spins up the star, enhancing the stellar 
wind, which draws more angular momentum from the 
system and increases the rate of orbital decay. The effect 
is that after a short amount of time the tidal spin-up 
fails to keep up with the angular momentum lost to the 
wind and the spin-orbit lock is lost. 

It should be noted that the above tidal evolution 
equations are only valid under the assumption of good 
alignment between the orbital and stellar angular 
momenta. This has been checked only for a subset of the 
currently known transiting extrasolar planets, through 
the measurement of the Rossiter-McLaughlin effect, 
and a non-trivi a l frac tion of misaligned pla nets is found 
(iJohnson et all 120 lit fWinn ct al. 201X1 : iNarita et al.l 



2010bllal: iWin n ct al.' 2010a,b; Bavhss ct al. 120101:^ 
Tripathi et alT 2010; Johnson ct al. 2009; Wi nn et al 



2009; 'Johnson et al. 2008; Winn et al. 2008; N arita et al 
200a i2007t lAVinn et al.i i2006t iQueloz et al.l 120001: 
Bundv fc Marcvl 120001: iSnellenI |2004 etc.). However, 



even for significant misalignment, the tidal torques will 
only change by a factor of order unity, much smaller than 
the uncertainties on Q*, which range over several orders 
of magnitude. In addition, the current measurements of 
the stellar spin-orbit alignment suggest that the shortest 
period planets, which are the only ones sensitive to the 
tidal dissipation, are well aligned with their parent star's 
rotation. 

Figure [1] shows an example of the evolution of HAT- 
P-20 b, computed with = 10^. The left boundary 
of all our plots has been placed at 30 Myr in order to 
show more details in the part of the evolution that is 
important for this analysis. 

From the top right plot one can see that for the first 55 
Myr, due to its pre-main sequence contraction, the star 
spins faster than the planet orbits. As a result the semi- 
major axis increases (to a degree not noticeable on the 
plot), and the star spins down, in spite of the fact that 
it is contracting. During the subsequent approximately 
35 Myr, the orbit and the convective zone of the star are 
tidally locked. This lock is quickly lost due to the star 
shrinking (bottom left panel) , causing the tidal coupling 
to sharply decrease. 

The core and the envelope of the star are clearly de- 
coupled for the first two Gyrs of the evolution. Initially, 



the core rotates slower than the convective zone, due 
to the fact that its moment of inertia does not change 
quite so much. However, at around 35 Myr, the situ- 
ation is reversed due to the stellar wind (and initially 
the planet) taking angular momentum away from the 
convective zone, but not the core. Significant differen- 
tial rotation is maintained to an age of about 2 Gyr by 
the stellar wind, which is gradually losing strength as 
the star spins down, until the core and the envelope are 
completely coupled. 

By around 3 Gyr, the orbit has shrunk enough for tidal 
torques to once again dominate the rotational evolution 
of the star causing it to spin up to a period of a few days 
as the planet inspirals. 

3. OBSERVATIONAL DATA 

From the more than one hundred currently known 
transiting planets, we based our analysis on fifty-three. 
Those were selected to orbit stars with masses between 
O.25M0 and 1.25M0. The lower cutoff was imposed be- 
cause we do not have reliable stellar models for masses 
below this range. For masses above the upper limit, 
the dissipation is likely dominated not by the convective 
zone, which at this point is next to non-existent, but by 
some dissipation mechanism in the radiative bulk of the 
star. Hence, assuming the same value applies beyond 
this point is not reasonable. 

In particular, probably the planetary system most of- 
ten given as an example of a very fast tidal orbital cvo- 
luti on, \yASP-18 b, (Hcllicr ct al. 2009; Hansen 201^; 
'Lai"20T2l; 'Hel lier et al.l i2011bl: iPenev fc Sasselovl l20lTt 
iBrown et al.ii2011l cf.) is not among the planets included 
in this work, due to the fact that its star lies above the 
1.25M0 cutoff we impose. The planetary orbit in this 
system, even just by itself, argues strongly against effi- 
cient tidal dissipation in the star. However, as discussed 
above, the mechanism of this dissipation is likely different 
than for the majority of the exoplanet host stars found 
by transit searches. 

Further, we restrict our sample to only systems which 
are consistent with having a circular orbit, and age lim- 
its quoted in the literature at least partially overlap with 
our 10 Gyr cutoff. Finally we exclude Kepler and CoRoT 
planets, because they are generally subject to much dif- 
ferent biases. Table [T] lists the systems and their relevant 
parameters that were included in our analysis, along with 
references to where those parameters were published. 

4. CORRECTING FOR OBSERVATIONAL BIASES 

For each transiting planet system, we need to figure 
out what the probability is that it would be observed 
and detected at any moment during its evolution. We 
split this probability into four parts: 

• pti'ansit. rpj_^g geometric probability that the planet's 
transit is observed from the Earth. 

• pdctcct. rpj_^g probability that the orbital phase cov- 
erage of a survey allows the detection of the transit. 
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followup . 



The probability that the transit candi- 
date will be chosen for follow-up by the survey and 
confirmed. 

• p'^se. 'pjjg distribution of ages of target stars for 
transiting surveys expressed as a probability. 
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Fig. 1. — Example of the evolution of the HAT-P-20 system for Q* = 10 . The various curves and points have the following meanings: 
a: the orbital semi- major axis in units of Rq; R,: the radius of the star in units of Rq; i?rad- the radius of the stellar core-envelope 
boundary in units of Rq; Porb- the orbital period in days; Pconv, Prad^ the spin periods in days of the convective and radiative zones of 
the star, respectively; Iconv, -'^rad- the convective and radiative zone moments of inertia in solar units; and p'r^^nsit^ pdctcct^ pfollowup^ page 
and Polder a^re defined in Section l4l 



The final probability density that we use is the 
normalized product of these four quantities. 

The transit probability (^transit -j ^^le simplest of the 
three biases. It is simply proportional to the ratio of the 
stellar radius to the orbital semi-major axis. 

The detection probability (pdetect-j jg more com- 

plicated, because the requirement for detecting a transit 
varies by system. For example, to detect a relatively 
deep event around a quiet bright star, observing only 
a few transits might be sufficient, while the detection 
of a shallower transit around an active faint star might 
require many transits (we incorporate the dependence 
of transit detection probability on stellar activity into 
pfoiiowup^_ In addition, p^ctcct jg ^ smooth function 

of orbital period, but rapidly oscillates, has sharp local 
maxima or minima near periods close to an integer 
multiple of 24 hr, e tc. fcf. |Collier Cameron et al.l 120061: 
ISmith et all 120061: IBurke et al.l 120061: IHartman et al.l 



1200 9b'). Instead of attempting to address all those 
complications, we will assume a simple smooth de- 
pendence of on the orbital period and present 
results with two different prescriptions for pd-<^*<^'^*_ The 
particular dependences of p'^'^*-'^'^* on orbital period we use 
are given in Figure [2] These roughly follow the curves 
published by various surveys for the recovery probability. 

The long period tail of p'^^*^<^^ is not well constrained, 
and is survey dependent (hence the two different 
prescriptions). However, it is also not particularly 
important, since planets with long periods are less 
affected by tides, even for relatively small values, 
and so our results are not sensitive to this assumption. 

The reason for prescribing shallower dependence of 
pdctoct pgj-jQ(j fQj. HAT survey is that, unlike all 
other transit surveys, HAT combines observations from 
two sites (one in Arizona and one in Hawaii) which 
increases their sensitivity at longer periods. 
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TABLE 1 

The Planetary Systems, their Parameters and the References used 



Planet Name 



I M, (Mo) R* {Rq) Stellar age (Gyr) 



Planet Mass Semi-major 



References 









Nominal 


Min 


Max 


V '-'J up J 


axis (AU) 


HAT-P-3 b 


0.917 


0.799 


1.6 


0.3 


4.5 


0.591 


0.03866 


HAT-P-5 b 


1.16 


1.167 


2.6 


0.8 


4.4 


1.06 


0.04075 


HAT-P-10 b 


0.83 


0.79 


7.9 


4.1 


11.7 


0.487 


0.0435 


HAT-P-12 b 


0.733 


0.701 


2.5 


0.5 


4.5 


0.211 


0.0384 


HAT-P-13 b 


1.22 


1.56 


5 
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0.851 


0.0426 
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0.0413 
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HAT-P-19 b 


0.842 
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8.8 


3.6 


14 


0.292 


0.0466 


WAT P 90 K 


n 7^(\ 

U. 1 0\) 




7 
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12 4 




U.UOUl 


HAT-P-22 b 


0.916 


1.04 


12.4 


9.8 


15 


2.147 


0.0414 


HAT-P-23 b 


1.13 


1.203 


4 


3 


5 


2.09 


0.0232 


HAT-P-24 b 


1.191 


1.317 


2.8 


2.2 


3.4 


0.685 


0.0465 


HAT-P-25 b 


1.01 


0.959 


3.2 


0.9 


5.5 


0.567 


0.0466 


HAT-P-26 b 


0.816 


0.788 


9 


4.1 


12 


0.059 


0.0479 


HAT-P-27 b 


0.945 


0.898 


4.4 


1.8 


8.2 


0.66 


0.0403 


HAT-P-28 


1.025 


1.103 


6.1 


4.2 


8.7 


0.626 


0.0434 


HAT-P-29 b 


1.207 


1.224 


2.2 


1.2 


3.2 


0.778 


0.0667 


HAT-P-32 b 


1.16 


1.219 


2.7 


1.9 


3.5 


0.86 


0.0343 


HD 189733 b 


0.82 


0.73 


7 


0.6 


7 


1.13 


0.03142 


OGLE-TR-56 b 


1.17 


1.32 


2.7 


2.6 


2.8 


1.29 


0.02386 


OGLE-TR-113 b 


0.78 


0.77 


0.7 


0.7 


10 


1.32 


0.0229 


OGLE-TR-182 b 


1.187 


1.53 


4.3 


2.4 


4.8 


1.06 


0.05205 


Qatar- 1 b 


0.85 


0.823 


6 


6 


13 


1.09 


0.02343 


TrES-2 


0.98 


1 


5.1 


2.4 


7.8 


1.253 


0.03556 


WASP-4 b 


0.92 


0.907 


5.5 


3.5 


8.7 


1.215 


0.02312 


WASP-5 b 


1 


1.084 


3 


1.7 


4.4 


1.637 


0.02729 


WASP- 10 b 


0.703 


0.775 


0.8 


0.6 


1 


3.07 


0.0371 


WASP- 13 b 


1.03 


1.34 


8.5 


3.6 


14 


0.46 


0.0527 


WASP-16 b 


1.022 


0.946 


2.3 


0.1 


8.1 


0.855 


0.0421 


WASP- 19 b 


0.97 


0.99 


5.5 


1 


14.5 


1.168 


0.01655 


WASP-21 b 


1.01 


1.06 


12 


7 


17 


0.3 


0.052 


WASP-22 b 


1.1 


1.13 


3 


2 


4 


0.56 


0.0468 


WASP-24 b 


1.184 


1.331 


3.8 


2.5 


5.1 


1.071 


0.03651 


WASP-25 b 


1 


0.92 


0.02 


0.01 


3.98 


0.58 


0.0473 


WASP-26 b 


1.12 


1.34 


6 


4 


8 


1.02 


0.04 


WASP-28 b 


1.08 


1.05 


5 


3 


8 


0.91 


0.0455 


WASP-34 b 


1.01 


0.93 


6.7 


2.2 


13.6 


0.59 


0.0524 


WASP-35 b 


1.07 


1.09 


5.01 


3.85 


6.17 


0.72 


0.04317 


WASP-37 b 


0.925 


1.003 


11 


7 


14 


1.8 


0.0446 


WASP-38 b 


1.203 


1.331 


5 


5 


14 


2.691 


0.07522 


WASP-39 b 


0.93 


0.895 


9 


5 


12 


0.28 


0.0486 


WASP-41 b 


0.95 


1.01 


1.8 


7 


7 


0.92 


0.04 


WASP-43 b 


0.58 


0.598 


7 


0.3 


7 


1.78 


0.0142 


WASP-44 b 


0.951 


0.927 


0.9 


0.3 


1.9 


0.889 


0.03473 


WASP-45 b 


0.909 


0.945 


1.4 


0.4 


3.4 


1.007 


0.04054 


WASP-46 b 


0.956 


0.917 


1.4 


0.8 


1.8 


2.101 


0.02448 


XO-2 


0.98 


0.97 


5.3 


4.3 


6.3 


0.57 


0.0369 


55 Cnc e 


0.905 


0.943 


10.2 


7.7 


12.7 


0.027 


0.0156 
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Hartman et al. (2009a) 
Bakos et al. (2009a) " 
Winn et al. (20l'Sa) 
Buchhave et al. ^2010) 
Hartman et al. (2011b) 
Hartman ct al. (' 2011b) 
Bakos et al. (2011) 
Bakos et al. (2011) 
Bakos et al. (2011) 
Kipping ct al. (2010) 
Quinn et al. (2012) 
Hartm ai^^taL ("2011a) 
Bekv et aH^ll) 
Buc hhave et al. (2011') 
Buchhave et al. (2011) 
Hartman et al. (2011c) 
Winn et al. (2006) 
Southworth (2010) 
Pont et al. (2007J 
Southworth (2010) 
Gillon et al. (2006) 
Southworth (201^ 
Alsubai et al. (2011) 
Sozzctti et al. (2004) 
Southworth (2010) 
Sozzetti et al. (2007) 
Daemgcn et al. (2009) 
Southworth (2010f^ 
Wilson et al. (2008) 
Southworth ct al. (2009a) 
Sanchis-Oicda et al. (2011) 
Anderson ct al. ( 2008) 
Southworth ct al. (20b9b) 
Christian et al. (2009) 
Skille n et al. (2009) 
Xister et al. (2009) 
Hellier et al. ( 2011b') 
Bouchv ct al. ( 2oTo) 
Maxtcd et al. (2010) 
Street et al. (2 010) 
Enoc h et al. (2011b) 
Smalley et al. f2oTo) 




In addition to orbital period the detection probability 
will depend on the brightness of the star and on the 
amplitude and frequency dependence of stellar vari- 
ability. However, since extrasolar planets are typically 
found around main sequence stars, the stellar luminosity 
does not vary much, and hence the stellar brightness 
dependence is mostly irrelevant for our purposes, and 
is ignored in our model. The dependence on stellar 



variability is generally complicated and difficult to 
quantify, so we include it as part of the foUowup prob- 
ability, and only in the general sense that high stellar 
activity is associated with high stellar spin frequency 
and hence the probability of a given transiting system 
being detected drops as the stellar rotation period drops. 

The follow-up probability (pf°ii°wup) ^-^^ ^^^^ 
cult to quantify, since it is subject to non-deterministic 
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4 6 8 10 12 14 
Period [days] 

default non-HAT modified non-HAT 

default HAT modified HAT 

Fig. 2. — Assumed dependence of pdotcct orbital period for 
one station surveys (non-HAT) and the two station HAT survey. 
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Fig. 3. — Assumed distribution of ages for the target stars of 
transit search projects. 



human evaluation and limitations due to follow-up 
resources specific to each project. In this work, we 
will idealize the situation and assume that it only 
depends on the rotational period of the star (Prot), or 
its projected equatorial rotation velocity w*sinz {v* is 
the the equatorial rotation velocity of the star and i is 
the angle between stellar rotational axis and the line 
of sight). Stars rotating fast have a smaller probability 
to be chosen for follow-up because the radial velocity 
precision will be limited by the rotational broadening 
of the spectral lines. In addition, a star's variability 
and hence a survey's ability to recognize a transit signal 
are correlated with the stellar rotation, so there is a 
bias against detecting transits for quickly rotating stars, 
which we did not include in pdctcct_ j^^^ |-j^jg work, we test 
the dependence of our results on pfoHow^P by assuming 
two different forms for it: (1) constant and (2) constant 
up to w*sinj < 20 km s^^, followed by an exponential 
decay with the follow-up probability reaching half its 
maximal value at v* sini — 40 km s^^. 

Finally, our presc ription for jp'^^'^ is bas ed on the dotted 
hne of Figure 3 of iTakeda et al.l ()2007D . In particular 
we use a piecewise linear approximation to their curve 
(see Figure [3]). In addition we also present results with 
page _ const. 

In the right panel of Figure [TJ we show the computed 
probability density functions pf^^nsit^ ^detect^ ^foiiowup ^^^^ 
page correspond to the evolution of for HAT-P-20 b 
presented in the other panels of the same figure. We also 
show the resulting PoUcr, computed as: 



-Polder (age) 



j-dcath ^transitpdctcctpfollowuppagc 



d(age) 



(8) 

In other words, Polder represents the probability to ob- 
serve a planet at its current age or older. As a planet's 
orbit decays, we have less and less chance to observe the 
planet before it is tidally disrupted. 

5. TESTING Q. VALUES 



In this section, we outline the procedure we use 
to determine if a given value of the Q* parameter is 
consistent with the observed set of exoplanet systems. 

We begin by finding an initial (5 Myr after stellar 
birth) semi-major axis for each planetary system, 
which after following the orbital evolution, according 
to Equation (P— [7]), to the present time, results in the 
observed value of the semi-major axis. 

We then continue the evolution until one of the follow- 
ing happens. 

I. The semi-major axis of the orbit falls below (1) the 
radius of the star or (2) the Roche radius for the 
tidal destruction of the planet, whichever comes 
first. If only tidal decay drives the orbital evolu- 
tion, we expect planets to spend very little time 
between their Roche limit and the stellar surface 
(when the former is larger), and so for our pur- 
poses the distinction between the stellar surface 
and Roche radius is unimportant: planets that 
have crossed their Roche radii are as good as gone. 



2. The star reaches the end of its main-sequence life- 
time. 

3. The system reaches an age of Tmax = 10 Gyr and 
neither of the above conditions has occurred. 

Having the complete time evolution of each system, 

we use pfansit^ pdetect^ ^foUowup^ g^^^ ^age f^.^^^^ Scction 

H] to calculate the probability that a random observa- 
tion throughout its lifetime will catch it at any given 
moment. Integrating this probability density from 
the present age of the system onward gives us the 
probability of observing this system no earlier than its 
present age - Poider- Hence, by definition, we expect 
Polder to be Uniformly distributed in the range (0,1). 
If we find it is not, this means that either we did not 
properly account for some observational or astrophysical 
bias or the orbital evolution computed is not correct. 
Assuming the various biases are handled correctly, and 
the orbital evolution equations (Equations ((I]) - 0) are 



°-8^ 02 04 06 08 1.0 

-folder 

— logi(,(5=6 ■ ■ ■ logioQ=8 uniform 

- - logioQ=7 - ■ logioQ=9 

Fig. 4. — Cumulative distribution function of Polder fof various 
values of Q*. 

appropriate, departures from uniformity will be due to 
a mismatch between the assumed Q* and the real one. 
So a KS test allows us to reject values of which are 
inconsistent with the currently observed population of 
extrasolar planets, around stars for which the dissipation 
is likely dominated by the convective zone. 

A plot of the calculated cumulative distribution func- 
tion of Folder (CDF(Foidcr)) appropriately corrected for 
observational biases and assuming the nominal values for 
the observed system parameters quoted in the literature 
for various values of is presented in Figure S) The KS 
test p-values corresponding to comparing those curves 
against a uniform distribution are shown as the red (+) 
symbols in Figure O 

There are several observational uncertainties that may 
affect our results. However, the uncertainty in the stel- 
lar ages (often a factor of several) dominates over the 
uncertainties of all other quantities. Thus, this is the 
only uncertainty we will account for. To demonstrate 
the possible impact on our conclusions, we compute cu- 
mulative distributions of Poidor with two assumptions for 
the actual ages of stars — the nominal ages and the lower 
end of the age range given by the corresponding publi- 
cation. The resulting p- values obtained by performing a 
KS test against a uniform distribution are presented in 
Figure O 

The shift to smaller values when the exoplane- 
tary systems are assumed systematically younger makes 
sense. A smaller present age of the system means that 
the interval between now and when one of the terminal 
conditions described at the beginning of this section oc- 
curs represents a larger fraction of the total lifetime of the 
system. To offset this, a smaller value of Q, is needed, 
shortening the future life of the system. 

Because the uncertainty is different for each system 
and because it depends on factors like the method for 
determining stellar ages, follow up instrumentation, stel- 
lar models used, etc., we cannot modify the expected 
distribution to include such an uncertainty. Instead one 
can show that if we prescribe a distribution for the age 
of each system, the statistical p- value corresponding to a 
given Q* is the expectation of the usual KS p- value over 



^° 6.0 6.5 7.0 7.5 8.0 8.5 9.0 

logio(Q*) 

+ nominal ages x minimal ages 



Fig. 5. — KS test p- values of the comparison of the cumulative 
distributions of Polder i calculated by assuming that all stars' ages 
are equal to the nominal values (red +) and the minimum allowed 
value (blue x). against a uniform distribution for various values of 
Q*. 



the distribution of ages. Since evaluating this expecta- 
tion requires taking an integral over a space which has as 
many dimensions as the observed transiting planets, we 
use a Monte-Carlo approach to calculate it. The precise 
procedure is as follows. 

1. Calculate the evolution of each system for a set of 
present ages covering the observationally allowed 
range. 

2. Draw a random age for each system from some pre- 
scribed distribution. 

3. Calculate the corresponding Poidcr and perform a 
KS test. 

4. Average the results of many such iterations to get 
the final p-value, which incorporates the uncer- 
tainty in the system ages. 

Since we do not know the appropriate distribution to 
assume for each system we will consider two options: a 
uniform distribution over the allowed interval of ages, 
and a normal distribution truncated at the age limits 
centered on the nominal age, with a standard deviation 
equal to one quarter of the given age range. 

6. RESULTS 

The procedure described in the previous section was 
performed six times with various assumptions for p^ctcct^ 
pfoiiowup g^^^ page discussed in Section IH for the two 
per system age distributions described above and finally 
we considered a case where we do not allow the convec- 
tive and radiative zones to rotate at different frequencies 
(i.e., Tc = 0). Table [2] lists the set of assumptions we 
consider, and the resulting KS test p-values are plotted 
in Figure m 

Evidently, only the assumption about how stellar ages 
are distributed over their observationally allowed ranges 
makes a noticeable difference. This is not altogether sur- 
prising considering that of all the assumed probability 
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TABLE 2 

The Various Assumptions used when Performing the K-S Tests for the Various Q, Values. 



Label 



^followup 



Age Distribution rc(Myr) 



Default 

Mod pdctoct 

Uniform pfoll°wup 
Uniform p^S'^ 
Normal age dist. 
Coupled 



Red (solid) curve in Fig. [5 
Blue (dashed) curve in Fig, 
Red (solid) curve in Fig. |2 
Red (solid) curve in Fig. 2 
Red (solid) curve in Fig. 2 
Red (solid) curve in Fig. |2 



Decreasing for v* s'mi > 20 km/s As plotted in Fig, 

2] Decreasing for v* sini > 20 km/s As plotted in Fig, 

Uniform As plotted in Fig, 

Decreasing for v* s'mi > 20 km/s As plotted in Fig, 

Decreasing for v* s'mi > 20 km/s As plotted in Fig, 

Decreasing for v* s'mi > 20 km/s As plotted in Fig, 



Uniform 
Uniform 
Uniform 
Uniform 
Normal 
Uniform 





10° 




10-^ 






alu 


10-* 


> 




a, 


10-^ 








10-= 




10-" 



6.0 6.5 



7.0 



8.0 8.5 



9.0 



+ default 

X mod pdetect 

o uniform pf°"°«'"P 



7.5 

logio(Q') 

V uniform p'"^'' 
□ normal age dist 
A coupled 



Fig. 6. — KS test p-values as a function of Q* for the range of 
assumptions detailed in Table [2] 



distributions this is the most poorly constrained obser- 
vationally. The core-envelope coupling timescale also 
has little effect, since allowing differential rotation or not 
never results in the star rotating faster than the planet 
except during a very short period during the pre-main 
sequence phase (like in Figure [T]) . 

All assumptions considered lead to very similar con- 
clusions: at the 1% level log^g Q* > 7. The different 
assumptions, fortunately, lead to differences only in the 
low p- value range of Q*. 

7. DISCUSSION 

The constraints on the stellar dissipation parameter 
we derive based on the transiting planet systems de- 
tected by ground based transit surveys (Q* > 10^) are 
inconsistent with constraints derived from observing the 
ci rcularization of b i nary s tellar systems in open clusters. 

iZahn fc Boucheti ()1989f) shows that, in order for tides 
to suppress the eccentricity in binary stars up to the ob- 
served circularization cutoff period, 10^ is required 
during the pre-main sequence phase, which corresponds 
to a p-value of 1% in our analysis. One can somewhat 
circumvent this marginal contradiction by assuming an 
evolu tion of Q^, as the ste l lar st ructure changes. How- 
ever, iMeibom fc Mathieul (|2005i ) find that in order to 
explain the observed rate of circularization during the 
main sequence phase even smaller values are required — 
Q, - 10^ 

One way to reconcile this apparent inconsistency 
is to note that the ratio of tidal frequency to stellar 
rotation frequency is very different for binary star 
circularization and tidal inspiral of a planet onto its 



star. For binary stars, the two components of the system 
are synchronized on a very short time scale (compared 
to circularization), so the tidal frequency is exactly 
twice the stellar rotational frequency. In the case of an 
exoplanet inspiral the tidal frequency is much higher 
than the stellar rotation. 

This different frequency might lead to two differ- 
ent me chanisms domi n ating the di ssipat ion i n the tw o 
cases. lOeilvie fc LinI (|20M . 120071) and (l2005allbl ') 
point out that if the tidal frequency is within a fac- 
tor of two of the stellar rotational frequency, inertial 
modes are resonantly excited in the star, which could 
lead to strongly enhanced shear and hence dissipation. 
Since the inertial mode frequencies are restricted to lie 
between —2a;* and 2a;*, this mechanism cannot oper- 
ate in the case of exoplanet systems. The currently 
favored mechanism of dissipation for low mass stars 
in the frequency r egime of ex o planetary tides is tur- 
bulent diss ipation CZahn' '19661; 'Zahn fc Boucheti I198S; 



Goldreich fc Nicholson 1977; Goldreich & Keeley 1977 
Goodm an fc Oh 1997; Penev et al.n2007l i2009b lai 12011 



Penev fc Sasselov 2011). This less efficient mechanism 



could result in dissipation efhciencies consistent with 
the constraints derived in this paper. In particular 
iPenev fc Sasselov! (|2011[ ). based on direct simulations of 
turbulent dissipation find Q* ~ few x 10^ to 10^, consis- 
tent with our results here. 

Finally, Schlaufman et all ()2010D argue that the sam- 
ple of Kepler planets favors 10^ < Q* < 10*", apparently 
outside our range. However, their analysis does not even 
consider Q* values above 10^ (other than infinity), and 
they do not derive the statistical significance of their lim- 
its, or the sensitivity of their result on the various as- 
sumptions included in their model (e.g., the conversion 
of mass to radius and the exoplanet population synthe- 
sis models they use). All this makes it difhcult to make 
firm statements about the (in)consistency of the two re- 
sults. The best way to address this would be to repeat 
our analysis for the sample of Kepler planets, but this is 
clearly outside the scope of this article. 

Our model suggests that the earliest stages of a close- 
in planet's dynamical history may be more complicated 
than widely considered. The top right panel of Figure 
[1] shows that, early on, HAT-P-20 rotated more quickly 
than its planet revolved and only after 35 Myr did the 
situation reverse. Consequently, the tidal torque exerted 
by the star switched signs at this point. If, for exam- 
ple, the planet were brought close-in through gas disk 
migration during its first tens of Myrs, then presumably 
there was a competition between torques from the gas 
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disk driving the planet in and tidal torq ues from the star 
driving the planet out. iLin et all ()1996l) pointed out the 
role that such tidal torques might play in stopping 51 
Peg b's inward gas disk migration but favored clearing 



of the gas disk very near the star for halting the inward 
migration. With so many more exoplanets in our sample 
now, many with much shorter orbital periods than 51 
Peg b, we plan to revisit this topic. 
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